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The effects of non-magnetic impurities randomly doped into a two-leg Heisenberg spin ladder 
are investigated. Using the continuous time quantum Monte Carlo loop algorithm we calculate 
the uniform and staggered susceptibilities of such a system. The obtained uniform susceptibility 
is well described in terms of an effective model of weakly interacting local moments induced by 
non-magnetic impurities for a 1% doping case, but not for higher concentrations. The staggered 
susceptibility however is significantly enhanced over that in the effective model already at 1% doping. 
Using a mean field approximation for the interladder coupling, we explain qualitatively the phase 
diagram of Sr(Cui_ I Zn a; )203. 



Spin ladder systems have attracted a lot of theoreti- 
cal and experimental interest recently, jj] The physics of 
these materials is believed to be closely related to that 
of planar cuprate superconductors. Theoretical studies 
have shown that the ground state of the undoped two- 
leg ladder is a short range singlet resonating valence bond 
(RVB) state with a finite spin-gap A k, 0.5J. This 
has been verified experimentally in materials such as 
SrCu2C>3 which has a gap of about 400K. ]2,|j Doping 
of holes into such a ladder system will probably lead to 
superconductivity of preexisting RVB singlet pairs with 
d-wave-type symmetries. This seems to be confirmed 
by the discovery of superconductivity in the ladder com- 
pound Sro.4Ca13.6Cu24O41.s4- 

In this letter we concentrate on another aspect, that of 
the suppression of the spin-gap by non-magnetic impuri- 
ties. This was first observed in the zinc-doped compound 
Sr(Cui_ x Z n;c ) 2 3 . |,§ Marked chan ges in the low tem- 
perature behavior occur upon partially substituting cop- 
per with magnetically inert zinc atoms. Susceptibility 
and specific heat measurements indicate that the impu- 
rities destroy the spin gap and lead to a long range Neel 
order at low temperatures. Doping with 1% impurities is 
enough to order the system magnetically at a Neel tem- 
perature of 3K. As the impurity density increases up to 
4% doping, the transition temperature increases to 8.2K. 
Further doping however reduces the Neel temperature. 

A theoretical model for this material is weakly cou- 
pled spin 1/2 Heisenberg ladders. Zinc doping leads to 
a random depletion of this ladder. Several aspects of 
an isolated depleted ladder have been studied. |7|-|l4| It 
was shown that a single impurity induces a local moment 
with S = 1/2 around it. f||l2[ The effective interaction 
between these two induced moments is ferromagnetic or 
antiferromagnetic depending on the sublattices of the two 
impurities, |l(],[l2j which does not contradict the Neel 
order. The magnitude of the interaction is determined 
by the distance between the two impurities, and decays 
exponentially. Based on this localize-moment picture, 



the enhanced staggered susceptibility has been discussed, 
jsjfhif but the actual magnitude of the staggered suscepti- 
bility, Neel temperature and its doping-dependence have 
not been established. 

Our main purpose in this paper is to calculate explic- 
itly the uniform and staggered susceptibilities and to dis- 
cuss the origin of the antiferromagnetic ordering and the 
phase diagram. Using the quantum Monte Carlo loop al- 
gorithm, fll5| , [l6[ we obtain the temperature dependence 
of susceptibilities, with high precision, down to extremely 
low temperatures. By analyzing these reliable numerical 
results, we can check whether the localized-moment pic- 
ture is valid or not for various values of impurity den- 
sity. It is found that such a picture is valid for only a 
1% doping case. Then, to determine the phase diagram, 
we consider the three-dimensional inter-ladder coupling. 
Actually, the interladder coupling inside the plane is frus- 
trated ferromagnetic El and its effect is not clear at 
present. Therefore in this paper we assume an unfrus- 
trated antiferromagnetic coupling which originates from 
the interplane coupling. This interaction is assumed to 
be weak compared to the intra-ladder couplings and we 
treat it in a mean field approximation. We estimate the 
staggered susceptibility of the three-dimensional coupled 
ladder system as 



Xmf 



1 - J'Xs 
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Here, \s is the staggered susceptibility of an isolated de- 
pleted ladder calculated by quantum Monte Carlo simu- 
lation and J' is an effective inter-ladder coupling. It is 
shown that, with this simple mean field approximation, 
we can explain the features of the phase diagram. In 
particular, we explain why a very small doping concen- 
tration is enough for the Neel order and why overdoping 
reduces the ordering. 

In this study we investigate ladders with randomly- 
depleted sites. So far, numerical calculations have been 
performed for periodically depleted spin ladders. [Hjll] 
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However we will show that the random distribution of 
impurities is essential for explaining the phase diagram. 
A typical configuration considered here is shown in Fig. 1. 
In the randomly-depleted case, the impurities can cut the 
ladder into disconnected non-interacting pieces as shown 
in the configuration in Fig. 1. Such configurations have 
not been studied before. In our quantum Monte Carlo 
simulations we have averaged more than twenty random 
configurations on a ladder of 1000 x 2 sites. In the fol- 
lowing discussion all the error bars include the statistical 
error of the Monte Carlo sampling and that arising from 
the different configurations of the disorder. Using this 
loop algorithm we can achieve Monte Carlo simulations 
down to extremely low temperatures such as T — 0.005 J, 
which cannot be achieved by other algorithms Thus we 
can discuss the temperature dependence precisely. 

First we discuss the uniform susceptibility. In Fig. 2, 
the temperature dependence of yT is shown for 1%, 5%, 
and 10% doping. Initially there is a Curie behavior of 
free spins at high temperatures (T ^> A), but we are not 
interested in this region and it is not shown in Fig. 2. At 
intermediate temperatures T < A and for a 1% doping 
case, we observe a plateau with \T = g 2 ii 2 B x/Ak^, with 
x being the impurity density. The horizontal lines in the 
figure are the lines with g 2 fj, B x j '4fce ■ This is consistent 
with the speculation that a local moment with S = 1/2 
is induced around an impurity. [^) In this case the mean 
distance between the two impurities is so large that the 
induced moments are essentially non-interacting in this 
temperature range. 

However at higher levels of doping such as 5%, we 
can see from Fig. 2 that such an independent localized- 
moment picture does not hold. In this case the interac- 
tion between the localized moments has already set in in 
this temperature range. 

Renormalization group arguments ]Io| for the random 
spin chain predict a limiting value of \T — ► g fj, B x/12kfe 
for T — > 0. These values are shown by arrows in Fig. 2. 
Apparently the predicted zero temperature limit is not 
observed for any impurity density down to T = 0.005 J. 
Since the localized moment picture seems to hold for 1% 
doping, the present result suggests that the renormal- 
ization group argument becomes valid at temperatures 
much lower than 0.005 J. On the other hand, for higher 
doping cases, the independent localized-moment picture 
is no longer valid, so that the temperature dependence is 
smeared out. In this sense, the renormalization group ar- 
gument does not seem to be relevant for the experiments. 
The results of the uniform susceptibility do not change 
within our statistical errors if we restrict the simulations 
to connected configurations. 

The staggered susceptibility (see Fig. 3) exhibits be- 
havior different from that of the uniform susceptibility. 
This is already visible in the undoped case, where the 
uniform susceptibility vanishes exponentially at low tem- 
peratures. The staggered susceptibility has a finite value 
at zero temperature. p7| 

Impurities lead to a divergent staggered susceptibil- 



ity at zero temperature. The magnitude of the stag- 
gered susceptibility at a high impurity density, depends 
strongly on whether we allow disconnected configura- 
tions or not. If we disallow such configurations we see 
a monotonous increase of the staggered susceptibility for 
a wide range of dopings. If we allow disconnected con- 
figurations we see that for more than 10% doping the 
staggered susceptibility decreases again. This explains 
the drop of the Neel temperature at higher doping levels. 
Impurities are needed to induce local moments. Over- 
doping however reduces these moments again by dilu- 
tion. We will discuss this effect more quantitatively later 
in this study. 

In Fig. 4 we plot the temperature dependence of XsT- 
One can see immediately that the picture of free spin-1/2 
moments induced by non-magnetic impurities does not 
describe the staggered susceptibility well even at 1% dop- 
ing. This behavior, quite different from that of the uni- 
form susceptibility, can be described as XsT = (a+b/T)T 
at low temperatures. The constant term a comes from 
the spins far from the impurities, and is almost the same 
as for the undoped case. The Curie term comes from 
spins around the impurities. However, the magnitude of 
this term is about 30 times as large as the corresponding 
term in the uniform susceptibility. The magnitude of this 
term at 1% doping can be explained by integrating the 
distribution of the excess z-component of spin around the 
impurity site, which is given by the DMRG technique, 
with a weakly staggered field. jlj] Thus, the antiferro- 
magnetic correlation between the impurities cannot be 
seen in our QMC method down to T/J = 0.005 at 1% 
doping. The crossover from the Curie behavior of in- 
duced spins to the coupled local moments is estimated to 
occur around 4% doping. 14 

Nagaosa et al. Q suggested that the staggered sus- 
ceptibility diverges as 
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(0.2) 



with a ~ 0.22, based on the renormalization group ar- 
guments of the induced-moment picture. However Fig. 4 
shows clearly that such a divergent behavior does not ex- 
ist at T — 0.005 J. This is for the same reason as the case 
of uniform susceptibility as discussed above: The renor- 
malization group argument is valid only in temperatures 
much lower than 0.005J. 

We argue here that, even if there is no strong diver- 
gence as in eq. (2), the significant enhancement of the 
staggered susceptibility leads to antiferromagnetic order- 
ing when the ladders are weakly coupled. Using the re- 
sults of staggered susc eptibility and the mean field ap- 
proximation of eq.(O.l) we calculate the Neel tempera- 
ture Tn from 



Xs(Tn) = 



(0.3) 



The exact magnitude of the inter-ladder coupling, J', is 
not known. We use the following estimates. The intra- 
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ladder coupling is taken to be J = 1000K as this gives 
the correct order of magnitude for the gap in the un- 
doped ladder. The inter-ladder coupling depends on the 



density as J' = J'°(l — x) 2 , where (1 



is the proba- 



bility of having two spins on adjacent sites. J'° is taken 
as 0.02J ~ 20K. This gives the same maximum Neel 
temperature as seen in the experiments and is consistent 
with the order of magnitude of the inter-plane coupling in 
other cuprate materials. In Fig. 5, we show the Neel tem- 
peratures as a function of impurity density estimated by 
this mean field approximation. As the impurity density 
increases up to 10% doping, the transition temperature 
increases to 6K. Further doping however reduces the Neel 
temperature. 

This phase diagram agrees qualitatively with the ex- 
periments on Sr(Cui_ x Zn :c )203 and can explain both the 
unexpectedly high Neel temperatures as well as the max- 
imum at substantially low levels of doping. The opti- 
mal doping ratio is however overestimated as compared 
to the experiments. We believe this to be an artifact 
of the mean field approximation which takes into ac- 
count neither the effect of disorder on the inter-ladder 
couplings nor the exact structure and frustration effects 
of the inter-ladder coupling in the plane. 

Note here that we can apply the mean-field approxi- 
mation for more than 5% doping where there is a cor- 
relation between the induced spins. For small impurity 
densities (less than 4%), the antiferromagnetic correla- 
tion between the impurities cannot be seen in our QMC 
method at temperatures as low as T/J = 0.005. There- 
fore the validity of applying the mean field approxima- 
tion for such small densities is not clear and the actual 
Neel temperature is expected to be much lower than the 
obtained value. In this sense, the question as to why 
high Neel temperatures are observed for small densities 
in Sr(Cui_ a; Zn a ;)203 is unanswered. 

In summary, we have found that randomly depleted 
spin ladders with weak three-dimensional interactions 
have antiferromagnetic ordering at low temperatures. 
The proposed speculation of independent free spin seems 
to be valid only at very low levels of doping of about 1% 
but fails at higher levels of doping. For higher impurity 
density, random depiction, which includes disconnected 
configurations in the ladder, is important to suppress the 
magnetic ordering. 

The authors thank M. Azuma for useful discussions 
regarding the experimental results of Sr(Cui_ x Zn x )203. 
The calculations were performed on the Hitachi SR2201 
massively parallel computers at the University of Tokyo 
using a parallel C++ Monte Carlo library developed by 
one of the authors (M.T.). 
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FIGURE CAPTIONS 

FIG. 1. Example of a randomly depleted ladder. 
Filled circles denote copper sites with S = 1/2 spin and 
open circles denote 5 = impurity sites. The configura- 
tion depicted here includes two impurities on neighboring 
spins which break the ladder into disconnected parts. 



FIG. 2. Curie constant per site for several impurity 
densities (1%, 5%, and 10%). The horizontal lines and 
arrows denote the predictions of the second plateau and 
the zero-temperature limit respectively for these doping 
levels, fllS 



FIG. 3. Staggered susceptibility per site as a func- 
tion of temperature. Results for the undoped system 
and for 1% and 5% doping are shown by triangles, open 
circles and closed circles, respectively. The inset shows 
results for 1%, 5%, 12.5%, 15% and 20% doping at low 
temperatures. 
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FIG. 4. Temperature dependence of % S T for 1%, 5% 
and 10% impurity doping. Note that the magnitude of 
the susceptibility for each impurity density is much larger 
than that of the uniform susceptibility. 



FIG. 5. Neel temperatures as a function of impurity 
density as estimated in the mean field approximation for 
the inter-ladder coupling. We have assumed J = 1000K. 



4 



o 



cut 

o 



o 



o 
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Fig. 2 Miyazaki et.al. 
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Fig. 3 Miyazaki et . al . 
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FIG. 5 Miyazaki et . al . 
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